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Abstract 

We present hadron-level predictions from the Monte Carlo generator Cascade and nu- 
merical calculations of charm and beauty production at the Fermilab Tevatron within the 
framework of the /c^-factorization QCD approach. Our consideration is based on the CCFM- 
evolved unintegrated gluon densities in a proton. The performed analysis covers the total 
and differential cross sections of open charm and beauty quarks, B and D mesons (or rather 
muons from their semileptonic decays) and the total and differential cross sections of bb 
di-jet hadroproduction. We study the theoretical uncertainties of our calculations and in- 
vestigate the effects coming from parton showers in initial and final states. Our predictions 
are compared with the recent experimental data taken by the DO and CDF collaborations. 
Special attention is put on the specific angular correlations between the final-state particles. 
We demonstrate that the final state parton shower plays a crucial role in the description of 
such observables. The decorrelated part of angular separations can be fully described, if the 
process gg* — > gg is included. 
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1 Introduction 



Charm and beauty production at high energies is subject of intense studies from both 
theoretical and experimental points of view p~H6]. From the theoretical point, these processes 
provide an opportunity to test the different predictions based on Quantum Chromodynamics 
(QCD) since the dominant production mechanism at high energies (i.e. small x) is believed 
to be quark pair production through the gluon-gluon fusion subprocess. At present, the 
problem of description of the charm and beauty production at the Tevatron within QCD 
is not fully solved. So, the difference between the DO and CDF measurements [IH1] of 
the 6-quark and 5-meson production cross sections and the calculations [7] performed in 
the framework of fixed-order next-to-leading logarithm scheme (FONLL) is about of a factor 
of 1.7. The central FONLL predictions [8] also lie below the data on the D , D + , D* + and 
cross sections which has been measured [5] by the CDF collaboration as a functions of their 
transverse momenta. Recently the calculations in the general-mass variable-flavour-number 
scheme (GM-VFNS) were performed [9j[l0] for the transverse momentum distribution of B 
and D mesons and have been found to be consistent with the data [HE]. However, there is 
still no GM-VFNS predictions for other measured quantities (like distributions in rapidity 
and azimuthal angle difference between the momenta of final mesons). 

Heavy flavour production has been considered [TTHT6] also in the framework of the kr- 
factorization QCD approach [17]. This approach is based on the Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) [18] or Ciafaloni-Catani-Fiorani-Marchesini (CCFM) [H] equations for the 
non-collinear gluon evolution in a proton and gives the possibility to take into account large 
logarithmic terms (proportional to lnl/x and, in the case of CCFM, also to In 1/(1 — x)). 
A detailed description and discussion of the /^-factorization approach can be found in [20]. 
A reasonable agreement between the /^-factorization predictions and the Tevatron data on 
the heavy flavour production has been found in [X2HT6] . It was demonstrated [TTlfT5lfT6] 
that studying the specific angular correlations between the transverse momenta of produced 
quarks can give an unique information about the non-collinear gluon evolution in a proton 
since taking into account the non-vanishing initial gluon transverse momentum in the 
/cT-factorization approach leads to the violation of back-to-back kinematics event at leading 
order. 

Measurements of bb di-jet total and differential cross sections (as a functions of the leading 
jet transverse energy and the di-jet invariant mass) and azimuthal angle correlations 
between two 6-jets have been performed by the CDF collaboration [21]. These measurements 
significantly extend the energy range investigated by previous analyses [1-5]. Studying the 
6-jet production is specially interesting since there are no additional assumptions on the 
fragmentation of the beauty quark into the B meson. The CDF collaboration has reported 
the preliminary data [6] on the charm pair production, where the D°, D*~ pair cross section 
and the D + , D*~ pair cross section as a function of the azimuthal angle between two charmed 
mesons have been meausred. The angular correlations in I?-meson production have been also 
measured [IH3] by the DO and CDF collaborations. 

The main goal of present paper is to give a systematic analysis of all available experi- 
mental data [TH(3ll2~T] on the heavy flavour production at the Tevatron in the framework of 
the /c^-factorization formalism. We produce the relevant numerical calculations in two ways. 
First, we will perform analytical parton-level calculations (labeled as LZ in the following) 
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similar to that done in [T51ITB] . In these calculations we will use the CCFM-evolved gluon 
densities (23] as default sets. The measured cross sections of heavy quark production will be 
compared with the predictions of Monte Carlo event generator Cascade [21] . Cascade is 
a full hadron level Monte Carlo event generator for ep, jp, pp and pp processes, which uses 
the CCFM evolution equation for the initial state cascade in a backward evolution approach 
supplemented with off-shell matrix elements for the hard scattering^]. In this way we will 
investigate the influence of parton showers in initial and final states for the description of 
the data. We will study the possible sources of theoretical uncertainties of our predictions 
(i.e. uncertainties connected with the gluon evolution scheme, heavy quark mass, hard scale 
of partonic subprocess and the heavy quark fragmentation functions). To investigate the 
dependence of our predictions on the non-collinear evolution scheme we will apply the unin- 
tegrated gluon densities derived from the usual (DGLAP-evolved) parton distributions (in 
the Kimber-Martin-Ryskin (KMR) [2E] approximation). Our special goal is to study specific 
kinematic properties of the final heavy quark-antiquark pair which are strongly related to 
the non-zero initial gluon transverse momentum. 

The outline of our paper is following. In Section 2 we recall shortly the basic formulas 
of the fc^-factorization approach with a brief review of calculation steps. In Section 3 we 
present the numerical results of our calculations and a discussion. Section 4 contains our 
conclusions. 



2 Theoretical framework 

The main formulas have been obtained previously in [T1T[T6] . Here we only recall some 
of them. The cross section of heavy quark hadropro duct ion at high energies in the kx- 
factorization approach is calculated as a convolution of the off-shell (i.e. ^-dependent) 
partonic cross section a and the unintegrated gluon distributions in a proton. It can be 
presented in the following form: 

cr(pp^QQX)= f — — -A(x 1 ,k 2 1T ,fi 2 )A(x 2 ,k 2 2T ,fi 2 )\M(g*g* ^ QQ)\ 2 x 

7 2 ji 2 ji 2 j > #1#2 ' 

x dp 1T dk 1T dk 2T dy 1 dy 2 — — , 

where A(x, k^, /i 2 ) is the unintegrated gluon distribution in a proton, \A4(g*g* —> QQ)\ 2 is 
the off-shell (i.e. depending on the initial gluon virtualities k\ T and k\ T ) matrix element 
squared and averaged over initial gluon polarizations and colors, and s is the total center-of- 
mass energy. The produced heavy quark Q and anti-quark Q have the transverse momenta 
Pit and p 2 T and the center-of-mass rapidities y\ and y 2 - The initial off-shell gluons have a 
fraction X\ and x 2 of the parent protons longitudinal momenta, non-zero transverse momenta 
ki T and k 2T (k 2 T = — k\ T ^ 0, k^ T = — k\ T ^ 0) and azimuthal angles <pi and <p 2 . The 
analytic expression for the \Ai(g*g* — > QQ)\ 2 can be found, for example, in [T^fT?]. 

The unintegrated gluon distribution in a proton A(x, k T , /i 2 ) in (1) can be obtained from 
the analytical or numerical solution of the BFKL or CCFM evolution equations. In the 
numerical calculations we have tested a few different sets. First of them (AO) was obtained 



1 A discussion of the phenomenological applications of Cascade can be found in [25] . 
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in [23] from the CCFM evolution equation. The initial (starting) distribution Aq(x, k^, /x 2 ,) 
has been used in the following form: 

xA {x,k 2 T ,fi 2 ) = Nx po {l - x) pi exp{-k 2 T /k 2 ), (2) 

where all input parameters have been fitted to describe the proton structure function F 2 (x, Q 2 ) 
An equally good fit can be obtained using different values for the soft cut and a different value 
for the width of the intrinsic kr distribution (the CCFM set BO). A reasonable description 
of the F 2 data can be achieved [23] by both these sets. 

To evaluate the unintegrated gluon densities in a proton A(x, k 2 -,, /x 2 ) we apply also the 
KMR approach [26J . The KMR approach is a formalism to construct the unintegrated parton 
(quark and gluon) distributions from the known conventional parton distributions xa(x, /J, 2 ), 
where a = g or a = q. In this scheme, the unintegrated gluon distribution is given by the 
expression 




where P a b(z) are the usual unregulated leading order DGLAP splitting functions, q(x, /x 2 ) and 
g(x, fj 2 ) are the conventional quark and gluon densities^ and T 9 (k^,/x 2 ) is the Sudakov form 
factor. The theta function 0(A — z) implies the angular-ordering constraint A = /x/(/x+ |kr|) 
specifically to the last evolution step to regulate the soft gluon singularities. 

In Fig.[Hwe plot the CCFM set AO (as a solid lines), the CCFM set BO (as a dashed lines) 
and KMR (as a dotted lines) unintegrated gluon densities at probing scale /i 2 = 100 GeV 2 
as a function of k\. for different values of x. One can see that the shapes of gluon densities 
under consideration are very different from each other In the following we will study the 
possible manifestations in the total and differential heavy flavour cross sections. 

In our analytical calculations (LZ) the multidimensional integrations in eq.(l) have been 
performed by the means of Monte Carlo technique, using the routine VEGAS [3U] • The full 
C++ code is available from the authors on requesto 



3 Numerical results 

We now are in a position to present our numerical results. First we describe our input 
and the kinematic conditions. After we fixed the unintegrated gluon distributions, the cross 
section (1) depends on the renormalization and factorization scales /Jr and \ip. In the 
numerical calculations we set /x^ = m| + (p 2 T + p| T )/2, fi 2 F = s + (where Qr is the 
transverse momentum of initial off-shell gluon pair), m c = 1.4 ±0.1 GeV, m;, = 4.75 ± 
0.25 GeV and use the LO formula for the coupling constant a s (/j, 2 ) with rif = 4 active quark 
flavours at A QCD = 200 MeV, such that a s {M 2 z ) = 0.1232. 

Numerically, we have used the standard GRV 94 (LO) [27], MSTW 2008 (LO) [28] (in LZ calculations) 
and MRST 99 [29 (Cascade) sets. 
3 lipatov ©theory, sinp.msu.ru 
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3.1 Inclusive charm and beauty production 

We begin the discussion by studying the role of non-zero gluon transverse momentum A>r 
in the off-shell matrix elements involved in (f ). In Fig. |2]we plot the differential cross section 
for cc and bb pair production as a function of transverse momentum p^, rapidity y® and the 
azimuthal angle difference between the transverse momenta of produced quarks Acj)®® at 
^/s = 1960 GeV. The solid histograms correspond to the results obtained according to the 
master formula in eq.(l). The dotted histograms are obtained by using the same formula but 
without virtualities of the incoming gluons in partonic amplitude and with the additional 
requirement kf 2 T < /i^. There are no cuts applied on the phase space of the produced 
quarks. As it was expected, in the back-to-back region A(jfl® ~ it both results coincide 
with each other. However, we find that a sizeable effect appears at low A(jfi® . Therefore 
the non-zero gluon transverse momentum in the hard matrix element is important for the 
description of the data at low and mediate A(ffi® . This effect is more significant for charm 
production due to smaller x. We have checked that the predictions between the LZ and 
Cascade calculations agree well at parton level. 

Now we turn to the transverse momentum distributions of charm and beauty production. 
In the case of inclusive 6-quark production, the transverse momentum distribution has been 
measured [1] by the DO collaboration and has been presented in form of integrated cross 
section (as a function of 6-quark minimal transverse momentum p^mm) at the total pp energy 
yfs = 1800 GeV for \y b \ < 1 (note that there is no cut on the rapidity of b). Our predictions 
are shown in Fig. [3] and compared to the data. We find good agreement between the LZ and 
Cascade predictions. We find a good description of the data when using the CCFM-evolved 
(AO) gluon distribution. The results obtained by using the KMR gluon density are rather 
close to the NLO pQCD ones [31] (not shown) but lie below the data. The difference between 
the CCFM and KMR predictions comes from the different behaviour of these gluon densities 
(see Fig. [1]) which is due to absence of small-x resummation in the KMR distributions. 

The CDF collaboration has measured [HE] the transverse momentum distributions of B + 
and several D mesons (namely, D°, D + , D* + and D+) with \y\ < 1 (where y is the meson 
rapidities in the center-of-mass frame) at yfs = 1960 GeV. Our predictions are shown in 
Figs. H]-[6]in comparison with the data [HE]. The fragmentation of the charm and beauty 
quarks into a B and D mesons is described with the Peterson fragmentation function [32] 
with €b = 0.006 and e c = 0.06. According to [33, 34J , the following branching fractions are 
used: /(& -> B+) = 0.424, f(c -> D°) = 0.582, /(c -> D+) = 0.268, f(c -> D* + ) = 0.229 
and f(c — > Df) = 0.084. The observed difference between the LZ and Cascade predictions 
is due to the missing parton shower effects in the LZ calculations. We address this point in 
more detail in Section 13.31 

Since the predicted meson transverse momentum distributions depend on the quark-to- 
hadron fragmentation, we have repeated our calculations for B + and D° mesons with the 
shifted values of the Peterson shape parameter e, namely % = 0.003 and e c = 0.03. These 
values are also often used in the NLO pQCD calculations. Additionally, we have applied the 
non-perturbative fragmentation functions which have been proposed in [7J[8j[35] an d which 
have been used in the FONLL calculations. The input parameters were determined [8"|I3"5] 
by a fit to LEP data. The results of our calculations are shown in Fig. [71 We find that 
the predicted cross sections (in the considered px region) are larger for smaller values of 
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parameter e or if the fragmentation function from [7JIB1E5] is used. However, the typical 
dependence of numerical predictions on the fragmentation scheme is much smaller than the 
dependence on the unintegrated gluon density. 

The DO experiment has measured muons originating from the semileptonic decays of 
6-quarks at y/s — 1800 GeV [lj for 4 < p^ < 25 GeV, \r]^\ < 0.8 (for both muons) and 
6 < m wt < 35 GeV, where r]^ is the muon pseudo-rapidity and m MAt is the invariant mass of 
the produced muon pair. In [2] the measurements are extended to the forward muon rapidity 
region, namely 2.4 < < 3.2. To produce muons from 6-quarks in the LZ calculations, we 
first convert 6-quarks into 5-mesons (using the Peterson fragmentation function with default 
value e& = 0.006) and then simulate their semileptonic decay according to the standard 
electroweak theory. The branching of b — > // as well as the cascade decay b — > c — > // are 
taken into account with the branching fraction taken from [34]. The predictions of the LZ 
and Cascade calculations are shown in Figs. M and M We find that our predictions with 
both CCFM-evolved unintegrated gluon densities describe the experimental data for both 
the transverse momentum and rapidity distributions of muons reasonably well. 

The calculated total cross sections of the 6-quarks, B and D mesons and their decay 
muons compared to the CDF experimental data [HtS] are listed in Table [TJ In Table [2] the 
systematic uncertainties of our calculations are summarized. To estimate the uncertainty 
coming from the renormalization scale /j,r, we used the CCFM set A0+ and AO— instead 
of the default density function AO in Cascade. These two sets represent a variation of the 
scale used in a s in the off-shell matrix element. The A0+ stands for a variation of 2 while 
set AO— reflects Hr/2. In all heavy flavour analyses studied here, we observe a deviation 
of roughly +10% for set A0+. The uncertainty coming from set AO— is generally smaller, 
but still positive. The dependence on the heavy flavour masses is investigated as well using 
Cascade. We varied our default values of m b = 4.75 GeV by ±0.25 GeV and m c = 1.4 GeV 
by ±0.1 GeV. All heavy flavour cross sections are observed to be lower, if the masses are 
enlarged and vice versa. 

We turn now to the investigation of the angular correlations between the produced parti- 
cles. As it was mentioned above, such correlations have been not studied yet in the framework 
of GM-VFNS scheme (which has been applied for transverse momentum distributions of the 
charm and beauty mesons only [9j[10]|). Experimental data on the azimuthal correlations 
in charm and beauty production come from both the CDF and DO collaborations. In the 
case of 6-quark production, CDF data [3] refer to the B, B azimuthal angle distribution 
measured at \y\ < 1, p T (B) > 14 GeV, p T (B) > 7.5 GeV and yfe = 1800 GeV. The DO 
data p] refer to the muon-muon correlation measured in the region of 4 < pj. < 25 GeV, 
\r]^\ < 0.8 and 6 < m MAt < 35 GeV at the same energy. In the case of charm production, the 
CDF collaboration has presented the measurement [6] of D°, D*~ and D + , D*~ pair cross 
section as a function of angle between the two charm mesons for the kinematic range \y\ < 1, 
5.5 < p T (D°) < 20 GeV, 5.5 < p T (D*~) < 20 GeV and 7 < p T (D+) < 20 GeV. Our predic- 
tions are shown in Figs. [10] and [11] in comparison with the data [HISIE]. We observe that the 
predicted shapes of azimuthal angle distributions are very different for different unintegrated 
gluon density functions. This is in a contrast to the cross sections as a function of transverse 
momenta or rapidities where all unintegrated gluon densities gave a similar behaviour. We 
can conclude that the cross section as a function of A<f) is very sensitive to the details of the 
non-collinear gluon evolution in a proton. As it was pointed out in [TT][T5l[T6] , such observ- 
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Source 


°\ B ) 


tjyjJ) 


a \ D ) 






CDF data [jib] 


2.78 ± 0.24 


13.3 ± 1.5 


4.3 ±0.7 


5.2 ±0.8 


0.75 ±0.23 


AO (LZ/Cascade) 


3.47/2.76 


13.34/9.31 


4.53/3.45 


3.81/3.19 


0.49/0.35 


BO (LZ/Cascade) 


2.54/2.02 


9.55/6.78 


3.27/2.54 


2.78/2.32 


0.35/0.26 


KMR (LZ/Cascade) 


1.33/1.16 


6.92/5.06 


2.40/1.92 


2.13/1.75 


0.29/0.22 


KMR (LZ, MSTW 2008) 


1.16 


6.04 


2.17 


1.79 


0.25 



Table 1: The charm and beauty total cross section in pp collisions at \fs = 1960 GeV. The 
measurement [H|5] was obtained in the central rapidity region (\y\ < 1) at 6 < pr < 25 GeV 
and 5.5 < pr < 20 GeV in the case of B and D mesons, respectively. 



Source 


a(B+) 




a(D+) 


a(D* + ) 




CCFM set AO 




2.76 


9.31 


3.45 


3.19 


0.35 


CCFM set A0± 




±10% 


±6% 


±6% 


±5% 


±9% 


CCFM set A0- 




±1% 


±4% 


±4% 


±3% 


±3% 


m h = 5.0 GeV, m c = 


1.5 GeV 


-8% 


-4% 


-4% 


-4% 


-3% 


rrib = 4.5 GeV, m c = 


1.3 GeV 


±11% 


±4% 


±1% 


±4% 


±11% 


e b = 0.003, e c = 0.03 




±7% 


±21% 


±23% 


±21% 


±25% 


Total 


4_16% 
^8% 


i_23% 
=t 4% 


i_24% 
£ 4% 


i_22% 
£ 4% 


i_29% 
^3% 



Table 2: Systematic uncertainties for charm and beauty total cross sections in pp collisions 
at yfs = 1960 GeV obtained with Cascade. The measurement [HE] was obtained in the 
central rapidity region (\y\ < 1) at 6 < pt < 25 GeV and 5.5 < pr < 20 GeV in the case of 
B and D mesons, respectively. 
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ables can serve as an important and crucial test discriminating the different approaches of 
the small- x physics. The CCFM-evolved gluon densities overshoot the data at A0 ~ ir and 
tends to underestimate them at A<p ~ 0. We observe that the peak at A</> — > is not at all 
reproduced by the CCFM unintegrated gluon densities; we address this point in more detail 
in Section [3T3l where we discuss the the process gg* — > gg. However the KMR gluon density 
has obviously a very different kx distribution (see Fig. [1]) and therefore provides a better 
description. 

3.2 bb di-jet production 

Recently, the experimental study of the bb di-jet production in pp collisions at y/s = 
1960 GeV has been presented [2~Tll2"2"] . The total cross section has been measured in the 
kinematical region defined by \yi\ < 1.2, \y 2 \ < 1.2, E 1T > 35 GeV and E 2 t > 32 GeV. The 
differential cross sections as a functions of the leading jet transverse energy and of the di-jet 
invariant mass have been measured with the azimuthal angle correlation between the two 
jets. The 6-jets are reconstructed with the JETCLU cone algorithm with radius R > 0.4. 
Our predictions are shown in Figs. fT2ldT4l We observe that at large Et and at large dijet 
invariant masses M the predictions fall below the measurement. Note, however, that in this 
kinematical region the quark induced subprocesses (such as qq —> bb) becomes important but 
are not taken into account here. At small and moderate Et and M, where gluon induced 
subprocesses play the leading role, the overall agreement of our predictions and the data 
is quite good. The measured A0 distribution has a significantly different tail at small A</> 
compared to the predictions; this will be addressed in more detail in the next section. 

3.3 Parton shower effects and additional gluon processes 

As shown above, the full set of the data on the heavy flavour production is in general 
reasonably well described. However the measured cross sections as a function of A0 show sig- 
nificant differences. These distributions are directly sensitive to additional parton radiation 
treated by parton showers or by other additional processes. In the following we investigate 
the influence of the details of the parton shower and the process gg* — > gg with subsequent 
branching of one of the outgoing (off-shell) g — > bb on the A0 distribution. For these studies 
we use the CCFM set AO. 

In Cascade, the initial state parton shower is angular ordered and based on the CCFM 
evolution equations. The final state parton shower is also angular ordered, but based on 
the DGLAP evolution equations. To investigate the dependence on the parton shower, we 
calculated the cross sections without parton shower, only initial state, only final state and 
both initial and final state parton shower, as shown in Fig. [15j Here, b dijet production 
is selected as an example. We observe a very small contribution of the initial state parton 
shower, since in fc^-factorisation the initial state parton shower does not influence the &t 
of the gluons (since it is determined from the unintegrated gluon density). However, a 
significant contribution comes from the final state parton shower. The prediction with full 
parton shower underestimates the back-to-back region of the azimuthal separation of the 
two b jets. This suggests that the final state parton showering generates too much gluon 
radiation, which causes less correlated jets as shown in Fig. [15] (a). To investigate this further 
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we changed the final state parton shower scale Q 2 max from Q 2 max = 4m^ to Q 2 max = m 2 . As 
shown in Figure [T5] (b), this leads to a higher correlation of the b jets and a very good 
description of the back-to-back region of the jets is achieved. 

At lower values of the azimuthal separation of the dijet system, higher order processes 
are expected to contribute significantly. Therefore, we repeated the b dijet analysis for the 
process gg* — » gg. In the matrix element calculation, which was performed in [35], one gluon 
in the initial state is on-shell and the other one is off-shell, as shown in Fig. [16j The heavy 
quark pair is then produced via parton showers in the final state as g — > bb. As shown in 
Fig. [13 (a), this process contributes significantly to the tail of the distribution very well, 
without any additional adjustment of parameters. Part of the contribution from gg* — > gg 
with g —7- bb is already included in the matrix element g*g* — > bb, so simply adding both 
contributions would result in double counting [37]. In [22J the influence of multi-parton 
interactions in the collinear frame is investigated, and it was found, that the contribution 
even in the small A$jj region is small, because a large Et of the jets is required. 

In conclusion, the azimuthal cross section measurements can be very well described by 
adjusting the scale for the final state parton shower to Q 2 max = m 2 and by including the 
process gg* — > gg with subsequent branching g — > bb, which in a collinear calculation would 
contribute only at NNLO. 

4 Conclusions 

We have studied charm and beauty production at Tevatron energies within the framework 
of the fc r -factorization. Our calculations are based on the CCFM-evolved unintegrated gluon 
densities in a proton. The analysis covers the total and differential cross sections of open 
charm and beauty quarks, B and D mesons (or rather muons from their semileptonic decays). 
The cross sections of bb di-jet production have been studied also. Special attention was put 
on the specific angular correlations between the final-state particles. Using full hadron-level 
Monte Carlo generator Cascade, we investigated the effects coming from the parton showers 
in initial and final states. Different sources of theoretical uncertainties have been specially 
studied. 

We obtain good agreement of our calculations for all observables and the recent exper- 
imental data taken by the DO and CDF collaborations. We have demonstrated, that the 
parton shower plays a significant role in the description of the cross section. We obtain a 
good description of the correlations once the higher order process gg* — > gg with subsequent 
g bb splitting is included. 
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Figure 1: The unintegrated gluon densities in a proton A(x, k 2 -, /i 2 ) as a function of at 
fj, 2 = 100 GeV 2 . The solid, dashed and dotted curves correspond to the CCFM set AO, 
CCFM set BO and KMR distributions, respectively. 
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Figure 2: Importance of non-zero transverse momentum of incoming gluons in charm and 
beauty production at the Tevatron. The solid histograms correspond to the results obtained 
according to the master formula (1). The dotted histograms are obtained by using the same 
formula but now we switch off the virtualities of both incoming gluons in partonic amplitude 
and apply an additional requirement 2T < \j? R . 
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Figure 3: The transverse momentum distributions of open 6-quark production at the 
Teavtron. The kinematical cuts applied are described in the text. The solid, dashed and 
dotted histograms correspond to the results obtained with the CCFM set AO, CCFM set 
B0 and KMR unintegrated gluon densities. The first column shows the LZ results while the 
second one depicts the Cascade predictions. The experimental data are from DO pQ. 
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Figure 4: The cross section as a function of transverse momentum of B + meson hadropro- 
duction. The kinematical cuts applied are described in the text. The left histogram shows 
the LZ numerical results while the right plot depicts the Cascade predictions. Notation of 
all histograms is the same as in Fig. [3j The experimental data are from CDF |3]. 
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Figure 5: The cross section as a function of transverse momentum of D° [(a), (b)] and D + 
[(c), (d)] meson hadroproduction. The kinematical cuts applied are described in the text. 
The first column shows the LZ numerical results while the second one depicts the Cascade 
predictions. Notation of all histograms is the same as in Fig. [31 The experimental data are 
from CDF 0. 
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Figure 6: The cross section as a function of the transverse momentum of D* [(a), (b)] and 
.D+ [(c), (d)] meson hadropro duct ion. The kinematical cuts applied are described in the text. 
The first column shows the LZ numerical results while the second one depicts the Cascade 
predictions. Notation of all histograms is the same as in Fig. [31 The experimental data are 
from CDF 0. 



17 




-* io 4 L 



Q h- 
Q. 

T3 



10' r 



10" 



- 1 1 1 1 1 1 ! 


! 1 1 1 1 1 1 1 


■ ■ 




L_j 
















1 



(b) 



Cascade 



> 10 J r 

o 



P io 2 L- 



10 - 



• CDF 

e h =0.006 

b 

. . . e u =0.003 



_l I I I I I I I I I I I I I I I I L_ 

10 15 20 25 



pB* [GeV] 



10 



15 



20 




Pt [GeV] 



15 20 

P?° [GeV] 



Figure 7: The dependence of our predictions on the fragmentation scheme applied. The 
soild, dashed and dotted histograms correspond to the results obtained using the Peterson 
fragmentation function with e& = 0.006 (e c = 0.06), e& = 0.003 (e c = 0.03) and the non- 
perturbative fragmentation functions from [7|8,35j. The first column shows the LZ numerical 
results while the second one depicts the Cascade predictions. Here we use CCFM set AO 
gluon density for illustration. The experimental data are from CDF [HE]. 
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Figure 8: The rapidity distributions of muons arising from the semileptonic decays of B- 
mesons. The first column shows the LZ numerical results while the second one depicts the 
Cascade predictions. The kinematical cuts applied are described in the text. Notation of 
all histograms is the same as in Fig. [3J The experimental data are from DO [2] . 
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Figure 9: The transverse momentum distributions of muons arising from the semileptonic 
decays of B-mesons. The left histogram shows the LZ numerical results while the right plot 
depicts the Cascade predictions. The kinematical cuts applied are described in the text. 
Notation of all histograms is the same as in Fig. [3j The experimental data are from DO [lj. 
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Figure 10: The distributions in azimuthal angle difference between the produced B- mesons 
[(a), (b)] and muons [(c), (d)], originating from their semileptonic decays. The first column 
shows the LZ numerical results while the second one depicts the Cascade predictions. The 
kinematical cuts applied are described in the text. Notation of all histograms is the same as 
in Fig. [31 The experimental data are from CDF [5] and DO [TJ . 
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Figure 11: The distributions in azimuthal angle difference between the produced D°, D*~ 
and D + , D*~ mesons. The first column shows the LZ numerical results while the second 
one depicts the Cascade predictions. The cuts applied are described in the text. Notation 
of all histograms is the same as in Fig. El The experimental data are from CDF [BJ. 
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Figure 12: The leading jet transverse energy distributions of bb di-jet production. The left 
histogram shows the LZ numerical results while the right plot depicts the Cascade predic- 
tions. The kinematical cuts applied are described in the text. Notation of all histograms is 
the same as in Fig. [3j The experimental data are from CDF [21J. 




Figure 13: The invariant mass distributions of bb di-jet production. The left histogram 
shows the LZ numerical results while the right plot depicts the Cascade predictions. The 
kinematical cuts applied are described in the text. Notation of all histograms is the same as 
in Fig. [31 The experimental data are from CDF |21j . 
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Figure 14: The distributions in azimuthal angle difference between the produced 6-jets. The 
left histogram shows the LZ numerical results while the right plot depicts the Cascade pre- 
dictions. The kinematical cuts applied are described in the text. Notation of all histograms 
is the same as in Fig. [3j The experimental data are from CDF [21] . 
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Figure 15: Parton shower effects for b dijet production in (a). The four lines represent 
full parton shower (solid line), no parton shower (dashed line), initial state parton shower 
(dashed dotted line) and final state parton shower (dotted line). The change of the parton 
shower scale from Q max = 4m^ to Q max = m2 for b dijet production is depicted in (b). In 
all cases the CCFM set AO is used and the data are taken from CDF [21] and CDF [6]. 
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Figure 16: The process gg* — > gg is shown schematically. The matrix elements are calculated 
in |36|. 
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Figure 17: The processes g*g* — > QQ (dashed line) and gg* — > gg for b dijet production in 
(a) and D°, D*~ meson production in (b). In all cases the CCFM set AO is used and the 
data are taken from CDF [21] and CDF [6J. 



25 



